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WAVE MULTIPLE SCATTERING BY A FINITE NUMBER OF UNCLOSED 

CIRCULAR CYLINDERS 


E.I. Veliyev and V.V. Veremey 


Introduction 


Study of the problem of wave multiple scattering by structures /3* 
which consist of a finite number of screens is of interest in the 
resolution of a number of practical questions in UHF electronics and 
antenna and measurement technology. 

The solution of this problem involves complex numerical com- 
puter calculations. 

A large number of studies [1-9 j have dealt with the problem of 
wave multiple scattering by a finite number of smooth scatterers. The 
problem is solved by both rigorous and approximate methods in these 
studies. The simplest approximate method of solution of the problem 
of multiple scattering by N smooth scatterers is to disregard the 
effect of the screens of the structure on each other. It is suffi- 
cient for the solution to knew the solution of the problem of multiple 
scattering by each of the scatterers separately . As experiments show 
however, this method leads to very significant errors even when the 
scatterers are at a considerable distance from each other. All the 
remaining approximate methods are reduced to more or less complete 
accounting for the reciprocal influence of the screens [7, 8, 23]* 

The problem of multiple scattering of a plane electromagnetic 
wave by a finite number of flat strips is discussed in [10]. When 
all the strips are in one plane, the problem is reduced to the solu- /k_ 
ticn of one integral Fredholm equation of the second kind. Smooth 
scatterers and flat tape do not have distinct resonance properties. 

•Numbers in the margin indicate pagination in the foreign text . "~™ 
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Study of the properties of structures formed of high quality screens 
capable of intensely scattering an electromagnetic field therefore 
arouses special interest . 

As was shown in [14], an infinite periodic array which consists 
of unclosed circular cylinders has qualitatively new resonance prop- 
erties which are of great practical interest. These resonances are 
observed upon excitation of the natural oscillations of the array by 
the incident wave [14], 

Any actual structure cannot consist of an infinite number of 
elements. Solution of the problem of wave multiple scattering by a 
finite number of unclosed screens in the form of cylinders with lon- 
gitudinal slits is therefore of special importance. 

We discuss two rigorous methods of solution of such problems. 
Following [3], we call the first rigorous method, which is conven- 
tionally used in study of the properties of N bodies [12], self con- 
sistent . The excitation of each cylinder by the incident wave and 
unknown fields of adjacent scatterers is analyzed here. This leads 
to the solution of systems of simultaneous equations or to systems 
of integral equations. Such systems usually are solved by computer, 
and the basic complexity is to obtain a system of equations which is 
convenient for numerical solution. 

We call the second method, first used by Schwarzchild [11] in 
studies of wave multiple scattering by slits, iteration [1]. This 
method is directly associated with the first one. This method es- 
sentially consists of successive study of the stages of wave scat- /5 

tering by the elements of the structure. It is assumed that the solu- 
tion of the problem of wave multiple scattering by a single element 
is known. The scattered field of the system is presented in the form 
of superposition of the scattered fields of the individual elements 
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and the scattered field of an arbitrary S-th element in the form of 
the sum of an infinite number of scattering orders 
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A scattered field of the first order is the scattered field of 
a single S-th scatterer. The scattered field of the second order 

is the response of the S-th scatterer to superposition of the 

scattered fields of the first order of the remaining scatterers, etc. 

It is easy to note that the boundary conditions of the total field 
are fulfilled on the surface of each scatterer in this formulation of 
the problem. The convergence of the series by scattering orders fol- 
lows from physical considerations and can be proved mathematically. 

These two rigorous methods, developed for structures made of 
smooth scatterers, can also be applied to study of resonance systems 
of a finite number of unclosed circular cylinders. The interaction 
between these screens is of special importance In study of the elec- 
trodynamic properties of such structures. The use of the iteration 
method makes it possible to estimate the extent of this Interaction. 

Other advantages of this method over the extremely frequently used 
[16] self consistent method also are discussed in the work. 

1. Solution of Problem of Wave Multiple Scattering by Finite Number 
of Scatterers by Successive Scattering Method 

A plane H polarized wave strikes a structure which consists of /6 
a finite number (N) of scatterers in the form of unclosed circular 
cylindrical screens at angle to the OY axis (see Fig. 1). 

The surface of the cylinders are Infinitely thin and ideally 
conducting, cylinder rauli are p s , the angular dimensions of the 
slits are 26 s and the orientation angles are £ s (S=l, 2, . . ., N). 

The generatrices of the cylinders are parallel to each other. It Is 
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required to determine the field scattered by the structure. 



Fig. 1. Total scattering cross section vs. wave dimensions of struc- 
ture, $ «0: 1. kd 1 =4.0 ka, 0 -0 *1°; 2. kd 2 =7-0 ka, 6 1 =9 2 =1° . 

We consider the case of nonoverlapping screens. We will find 
the solution of the problem by the successive scattering method pro- 
posed in [1, 11], which has a number of advantages because of the 
clear physical interpretation. 


We write the total field of the structure in the form of super- 
position of the incident field and the scattered fields sought from 
each unclosed cylinder: 




W/ oa i Hi * Z?H, 


scat 




( 1 ) 


WQ.VG 

Field H everywhere outside the elements of the structure satis 



fies the Helmholtz equation and, on the surfaces of the cylinders, 
the boundary condition of Neumann, the condition of radiation to 
infinity and the condition of finite energy in any bounded region of 
space. 


Following [1], we present the scattered field of the S-th cyl- 
inder in the form of an infinite sum over the so called orders of 
scattering: 


u scat - 


E" 


( 2 ) 


We define the field of the first order of scattering as the response 
of a single S-th cylindrical screen to an incident plane wave, i.e., 
the Helmholtz equation, the radiation condition and the Neumann bound- 
ary condition of the surface of the S-th cylinder are satisfied by 
superposition of fields H ^+^H The scattered field of the first 

Zi Z 

order can consequently be written in the form [17] 



( 3 ) 


where (R c , 40 are polar coordinates centered on the axis of the 

^ b (1) S 1 

S-th cylinder, J' , H^ are Bessel functions, and are the 

Fourier coefficients of the current density function of the first 

order on the surface of the S-th cylinder. 


S 1 

For determination of coefficients \i , we subordinate total 
i S 1 

field H + H to the Neumann boundary condition on the surface of 
z z 

the S-th cylinder 



( 3 f ) 


In the system of functional equations of the first order obtained in 
this manner, inversion of the static portion of the matrix can be 
carried out by the method of the Riemann-Gilbert problem, and a sys- 
tem of linear algebraic equations of the second kind can be obtained 
[ 19 ]: 

ORIGINAL PAG£ 13 
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where 


M* " * *6* ' 


r i fs t **irfxfyfrietrt- c, «.««« 


C'- 

ntn 


n+O 


[ " ■ Wo , tn-O,n-0 ; 


tpi 


4 •> P 1 *) 

-<*M/ 4j^T ,*-< 

/t£= e 1 **}* Sin (Vi*y& 


n*o 


W 


(5) 


( 6 ) 


J is a random number from 1 to N. /9 

System of linear algebraic equations of the second kind (4) is 
Fredholmian [19] and can be solved by the reduction method. Its solu- 
tion can consequently be obtained with any preassigned accuracy. Sys- 
tem (4) can be solved by the successive approximation method In a 
number of cases [17], which substantially simplifies the calculations. 

By summing the scattered first order fields of each cylinder, we 
obtain the solution of the multiple scattering problem in the single 
scatterer approximation. As has been stated, this approximation gives 
a strongly distorted concept of the field scattered by the system. 

The degree of accuracy of the results can be estimated by means of 
accounting for the subsequent scattering orders. 

We will call the response of the S-th cylinder to superposition 
of the first order scattering field the second order scattering field 
of this cylinder, i.e., we assume that superposition of the scattered 
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first o^der fields from the cyl^der^wyt 1 ^ numbers i^S (1*1, 2, . . . , N) 
and the scattered second order field of the S-th cylinder satisfies the 
Helmholtz equation, the boundary conditions on the surface of the cyl- 
inder and the radiation condition. Vie write this total field in the 
form 


XV-U' ■%' 


q p 

where H is the second order scattering field of the S-th cylinder, /10 
z 

and the £’ symbol specifies that the value of summation notation q=S 
is excluded. 

The second order field which satisfies the wave equation and the 
radiation condition outside cylinder S (like the first order field) 
can be written in the form 

%•(*.») -Lpi Z.(*A)HS(**s)er* ( 7 ) 

(m) 

S n2 

With Eq. (3) and (7) taken into account, we write field H z 

(8) 

• flul X (Kfid ■HMe.ye™* 


S 2 S n2 

For determination of p , we subordinate field to the 

Neumann boundary condition on the surface of the S-th cylinder, i.e.. 


dg * l*i 'A 

*/jl* - o , * 


In subordination of the total field to the boundary conditions, 
a change must be made to the (R g , $ g ) coordinate system connected 
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with the S-th cylinder. The summation theorem for Bessel functions 
must be used for this 

( 10 ) 

Series (10) converges with R 0 <r„ o , <fc is the characteristic position 

o QS CJS 

angle of the q-th cylinder relative to the S-th cylinder and r qg is 
the distance between these cylinders. 

By using Eq. (10), we obtain from Eq. (9) the system of func- 

S 2 

tional equations for determination of y^ 


(11) 

(m) ' 

; eim< * • ° . *&**.<*. , 

where *Q‘ m , /II 


The Fourier coefficients of the first order surface current 
density function q p m 1 are determined from the solution of system 
( h ) . System of Eq. (11) is similar to the system of functional 
equations obtained by solution of the problem of excitation of a 
single cylinder by a plane wave, i.e., by determination of the first 
order scattering field. These systems only differ on the right 
sides, which characterize the exciting field. This is completely 
understandable, since both the first and second order scattering 
fields are defined as the response of the same cylinder to different 
perturbing fields. This is the incident plane wave in the first case 
and the first order scattering field of the remaining cylinders in 
the second. 

By applying the method of the Riemann-Gilbert problem to system 
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(II), we obtain a system of linear algebraic equations of the second 

3 p 

kind for y m , the matrix of which equals the matrix of system of 
Eq. (4) 



*Z c& W * V - 

(<*) ' 

e" (, ‘ v ;: , 

n 

-trfLZl 


( 12 ) 


Where . Ol 

System (12) is analagous to system (4). It can also be solved 
by the reduction method and by the successive approximations method 
in particular cases. 

We will call the response of the S-tb cylinder to the second 
order scattered field of the remaining cylinders the third order 
scattered field 

£'%' '‘Hi ■%" 

f* 

in which With ?s* 2* 

The third order scattered field is found just like the second 
order field was found. 


In a similar manner, it is easy to determine the fields of all 
the remaining orders of scattering. By then summing the scattering 
orders, we obtain the fieDd scattered by the structure 


*•/ fi (ml * 
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Q — 

°P m 1 are determined for random parameters of the structure by solution 
of the system of linear algebraic equations of the second kind 


where 


f e^'\7 , — 

Li ,uo. 


'q: 


d^o 


e"”* ,£.t 


* 7 ZZH *#. 




The resulting solution of the multiple scattering problem can /13 
be called rigorous if, by assigning the order of iterations, i.e., 
by taking account of a specific number of scattering orders, the solu- 
tion of the problem can be obtained with any preassigned accuracy. 

We note that the matrix of the system of equations is defined 
only by the cylinder number and does not depend on the scattering 
order. This substantially speeds up the calculation process and 
simplifies its management. 

More than that, in the particular case when the cylinders have 
narrow slits (6 s -0, S«l, ...» N), system (14) permits analytical 
solution. Analysis of the behavior of the coefficients of the matrix 
or the system of equations shows that, with Increase of indices n and 

c 

m and when 0g (angular dimension of the slit )-*■(), C nm with m^O, n^O, 
mj*n are very small, and they can be disregarded in solution of the 
system of equations, i.e., 

y. - cty -‘c, 

/n+0 , * Lmm • 
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In this manner, by leaving the diagonal elements and the elements 
of the zero column and zero line in the matrix of the system of equa- 
tions, the solution can be written in the form 

r r c*. ‘c/0-.z.) *‘t: 

< .*• *-£*<&*» (15) 

'/Urn * {jUf eZ. -Kyo-Cil , m*o. 


A numerical test of Eq. (15) was carried out. The results were /Ik 
compared with the rigorous solution of the problem of wave multiple 
scattering by a single unclosed cylinder. With e<10°, the error with 
ka<5*0 is from 0.01? to 1?. At angles 0<1° , Eq. (15) gives an ac- 
curate result (to within 1 ?) over a significantly wider range of 
change of ka (see Table 1). 


Because the solution of the multiple scattering problem as 60 +O 
of an unclosed cylinder Is written in explicit form, a closed expres- 
sion can be obtained for the field scattered by the structure by in- 
troducing some designations (similar equations were obtained in [ 1 ], 
which dealt with study of the scattering properties of circular cyl- 
inders ) 


Hr 


■ft &&) »?(**•) 

*•/ *»••• 




/ 


where 


C*'*) 


g p* • fl 21 21 3m* (y*» 0 ) (%****• Kf * 0 S 0 *) 


s 0 m* 


■A a ., 


•Zf 


e M r, 


in which 
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where 


X- 


V." . 

- l * (*fif > m ‘° 


J - C, 


i” 




tn + O 

CJL> Cm<9 


m'?9 * t-m'm' 


HI ~0 


/n 


l.i order to better explain the iteration procedure of obtaining 
the .solution of the problem of multiple scattering by a finite number 
of scatterers, we repeat the stages of calculation in order t6 show 
that the boundary conditions are fulfilled on the surface of each un- 
closed cylinder. 


We show that 



with £ s *Qf < (we will subsequently imply this condition). 

,c*j, 


In this manner 


A 4 Hi h ?~ 
j§g t H* H?J+Z! i d* s ' Wlth 




That is, the total field consists of the superposition of the inci- 
dent wave, the wave scattered by the S-th cylinder and the waves 
scattered by the remaining cylinders. By presenting the scattered 
fields in the form of an infinite sum of the scattering orders, we 
write 


d. ...tot /£ 
d*' Hl '(a?. 




with 
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TABLE 1. 
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k 
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1 
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0.3745 
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0.0418 

-0.0528 
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-0.1012 

-0.0286 

gg j 

1“ 
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0.3745 

0.1849 

0.0418 

-0.052b 


B 

1.5938 

-0.II99 

-o.ior 
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§■!■■■ 
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\ 

* 
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« 

$ 

V 
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S~ 
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■■■ 

mamm 
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so 

1.0493 

0.4876 

0.3796 

0.0508 
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-0.2515 

-1.5575 
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■0.3924 

-■% ■ 
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-1,5630 

1.0497 
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<r.tetrr ' 

Tm 
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1 

(a 
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Tm 
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■■ 

■■■ 


■■■ 

■BBBflHI 

■■■■■ 

mammal 

• 

$ 

<1 

<*> 

rv s 

p a 

Me 

-0.3221 

0.2313 

0.2045 

0.0131 

-0.0684 

c 


-1.0646 

-0.7286 

-0.3581 

-0.1633 

-0.0166 

V 

s* 

b 

i 

ib 

Me 

-0.3218 

0.2293 

0.2037 

' 0.1298 

-0.0686 

£ 

1 

Jm 

-1.0659 

-0.7207 

-0.3527 

-0.1617 

0.0156 



■■ 

■HI 

■Mi 

■■■Ml 

■■■■1 


1 

? a 

/ft- 

-1.6206 

-0.5715 

1.5045 

0.2148 

-0.3552 

1 

Tm 

-0.6306 

-4.7016 

-1.3099 

-1.4946 

-0.2277 

I 

5 b 

Me 

-I. 6601 

-0.8767 

1.4465 

0.1380 

-0.3180 

I 

Tm 

-0.6244 

-4.4720 

-1.4337 

-1.4337 

-0.1891 


Key: a. Rigorous b. Approximate 






















































































































ORIGINAL PAGE 19 
OF POOR QUALITY 

Because of equality (3 ? )» the sum of the first two terms equals zero /l8 


J 


with • * 7 *. 


9 jr'y^u* 
0*s4*£t * ' 


It follows from equality (9) that we can drop out the first two 
terms and in the same manner for all scattering orders, since the 
scattered fields from the remaining cylinders constitute the excita- 
tion of any individual cylinder and as a response to which the field 
of the next scattering order is formed. 


Consequently, the n-th order of excitation 


and hence 


d 





■£-Urj ■ O (in metal). 

d *s r Its -J>, 


The physical principles which are the basis of the Iteration al- 
gorithm of construction of the solution of the problem of multiple 
scattering by a finite number of scatterers give the right to speak 
of the convergence of the iteration process, i.e., the series by 
scattering orders converges This also can be proved 

-»w 

mathematically . 

It Is completely evident that the convergence rate of the 
series depends on the incident wavelength as well as on the geometric 
dimensions of the structure (primarily on the distance between scat- 
terers ). The rate of convergence of this series can be estimated nu- 
merically (see Section III). Estimates show that, even when vhe 
cylinders are located quite close to each other (there is one diameter 
between cylinders), the scattering order series converges quite ra- 
pidly (eight members of the series is sufficient to obtain results 
which are accurate to within 0.1$). 

The availability of clearly defined stages of construction of 111 
the solution of the problem and the algorithmic nature of the solu- 
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tion makes the iteration method very convenient for computer program- 
ming of the problems. Compared with direct self consistent method 
(described in Section II), there are substantial savings in working 
storage here, which makes rigorous solution of such problems possible 
even on a computer of medium capacity. A great advantage of the 
iteration method is the physical foundations of presentation of the 
scattered field in the form of superposition of the scattering orders. 
This method makes it possible to estimate the extent of interaction 
between structural elements and to show the region to which a given 
approximation can be restricted and the frequency region and struc- 
tures for which single scattering theory can be used, i.e., to gen- 
erally disregard the interaction. 


We note in concluding this section that the field within the 
cylindrical screens can be determined in the following manner: 


'tt-i 


where 


II. Rigorous Solution of Multiple Scattering Problem Based on Self 
Consistent Method 


Following [3], we will call the method described in this section 
self consistent. The substance of the self consistent method is that 
the field scattered by some cylinder is found on the assumption that /20 

the field scattered by the remaining cylinders is known, i.e., by 
considering all the scatterers In turn and presenting the exciting 
field in the form of the superposition of a plane wave and unknown 
scattered fields of the remaining elements of the structure, after 
application of the boundary conditions, we arrive at a system of 
simultaneous equations for the unknown amplitudes of the scattered 
fields of the structural elements [12]. 


Quite high efficiency of the self consistent method is ac- 
hieved in the case in question by partial Inversion of the matrix 


15 
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of the system of functional equations of the first kind, similar to 
the way this was done In study of the unclosed cylinder array [14], 

We present the total field in the form of superposition of the 
incident and scattered fields. The field scattered by the structure 
is written in the form of superposition of the fields scattered by 
the individual elements 


* / 4 sc ^//; *£% s ? at 


(16) 


We will use the (R g , $ g ) coordinate system (S=l, 2, . . . , N) tied 
to the centers of the corresponding cylinders to record the incident 
and scattered fields. 


By presenting the scattered field In the form of superposition 
of the bilayer potentials and by using the property of periodicity 
of the surface current density function on the cylinders with respect 
to <j> s , we obtain [17] 


nr-lZ &( wK'Kj , «, >a. 

fi (*; ' 


By applying the summation theorem for cylindrical functions, the 
field scattered by the structure can be written in the (R g , 4> s ) co- 
ordinate system of a random cylinder of number S 


..scat « . 


f- 



*<P rt _ 




( 17 ) 
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Total field (16) should satisfy the wave equation, radiation condi- /21 
tion, the condition of finite energy in any finite space, as well as 
the boundary conditions on the surface of each scatterer. 


The Neumann boundary condition applies in our case: 


A u t0 *l . 


with 

Sm 1,2, ..Ji. 


(18) 


The following condition is fulfilled in the remaining portion of the 
interval of change of <$> s : 

fitf e*** *0 , £-%<%<&'%• < 1 9) 

*—• 

We substitute the expressions for the incident and scattered fields 

(16), (17) in Eq. (18) and we differentiate. The resulting equation 

here together with Eq. (19) form a system of functional summator 

equations with a trigonometric kernel for the unknown y ^ : 

P 


where 



t* g 

e **Q» Si *(n *%*) 


( 20 ) 


/22 


A similar system of equations was obtained in [ 1 4 ] in solution of 
the problem of multiple scattering of a plane wave by an Infinite 
periodic array of unclosed circular cylinders. However, in view of 
the fact that the surface current densities on the cylinders are 
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different in the problem in question in distinction from [15, 14], 
the system of functional equations for y n q , with q=l, 2, . . N, 
must be considered Jointly. As is well known, systems of the first 
kind are unsuitable for numerical analysis. We therefore partially 
invert the matrix which corresponds to the problem of wave multiple 
diffraction by a single cylinder by the Riemann-Gilbert method. The 
procedure of expansion of the matrix in the completely continuous and 
easily inverted portion is discussed in detail in [20]. 


Thus, by using the Riemann-Gilbert method, we obtain a set of 
simultaneous linear algebraic equations of the second kind which are /23 
convenient for computer utilization: 





( 21 ) 


S S SSI 

where C and b were determined in Section I (with b = b ) 

nm n n n 


<• 


«o 

-&(*p*)* zU*j>s) ‘ 


The scattered field should satisfy the condition on the edge, i.e., 
the condition of limited energy in any finite volume of space. It 
can be shown that, to satisfy this condition, the Fourier coefficients 
of the surface current density function must be found in Gilbert 
space Jl 2 , i.e., the inequality 

ILwm 

should be fulfilled. 
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By using estimates for Bessel functions of higher order, it can be 
shown that 


£ WH 


< » 





System of Eq. (21) can be solved by the reduction method with /24 
any preassigned accuracy. 

The matrix of the resulting system of equations is of a distinct 
cellular nature. In some particular case of mutual location of the 
cylinders, the matrix of the system of equations can be cellularly 
Teplitsian, which permits the use of the efficient numerical algo- 
rithms developed for the calculation of such systems [18]. 

We discuss the total scattering cross section of a system of 
unclosed circular cylinders. 


We write the scattered field in the form 


H, Sca -2T //' , 


m £ jii tCfitV-e 4 "* 

~ / wi th /e, >r„ y,' 

A ’A %(«/>,) ^r p )e^ 


where 


By determining the flow of the Umov-Poynting vector through a cylin- 
drical surface of infinitely large radius with the center on the axis 
of one cylinder, we obtain 

0 <:■ eeje r H; *. </?. • ££/#•/'■ ( 22 ) 

o 

% 3 

The series found by calculation of y m proves to be slowly converging 
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however. Thus, Eq. (22) proves to be unsuitable for the calculations. 

By using optical theorem, a convenient equation for calculation 
of the total scattering cross section was successfully obtained: 

*'• *(*-**/!:)• 

Study of the total scattering cross section permits analysis of the 
scattering power of the structure. 

III. Discussion of Calculation Results 


Computer calculation of various characteristics of a structure 
which consists of a finite number of unclosed cylinders was performed 
by both the self consistent and iteration methods. The advantages of 
the iteration method was shown in the preceding sections. They are 
primarily the algorithmic and graphic nature of the method and then 
the economical use of the working memory of the machine, which makes 
it possible to calculate multielement structures on medium capacity 
machines without the use of external memory, and finally the advantage 
in calculation time. 

The question arises however of the convergence rate of the itera- 
tion process. It is physically obvious that the convergence is the 
higher the greater the distance between scatterers. 

Numerical computer utilization of the iteration method makes it 
possible for us to estimate the convergence rate of the successive 
scattering process. 

We limited ourselves to study of structures which consist of two 
scatterers in the calculations, since the purpose was to test the suc- 
cessive scattering procedure for subsequent application to study of 
multielement structures. 
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The calculation results were compared with the results of 
rigorous calculation of a similar structure by the self consistent /26 
method. It turned out that the total scattering cross section of a 
structure with kd^-4.0 ka and kd 2 =7.0 ka (where d is the distance be- 
tween scatterers, a is the cylinder radius, k-2 v/\) with eight suc- 
cessive scatterings taken into account can be calculated to within 
0.5/S (see Fig. 1-3 and Table 2). As should have been expected, the 
accuracy of calculations with kd^ is higher than when kd 2 , and 5-6 
iterations are sufficient for calculation of such a structure. The 
accuracy of calculations in the resonance regions is reduced because 
of amplification of the effect of the elements on each other. The 
calculation error is less than 0.5% however. 

The real and imaginary parts of the Fourier coefficients of the 
surface current density in response to the electromagnetic field next 
to the scatterers are presented in Tables 3-8. The Tables 3-8 data 
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are accurate to within 10 ’^. It is easily seen here how, by incor- 
porating the additions which respond in turn to all stages of scat- 
tering, i.e., repeated successively scattered waves, we arrive at a 
more complete accounting of the interaction between the scatterers. 



Fig. 3 . Total scattering cross section vs. ka: 1 . kd» 

2.1 ka, £,*£0*90°, 6 ,-e p= 5 °; 2. kd«2.1 ka, £,«0°, 

180 0 , e^e^ 0 . 

Thus, to obtain Fourier coefficients with less than 1 % error and 
ka- 0.4 and ic^, it is sufficient to allow for six successive scatter- 
ings and, for ka» 0.4 and kd 2 , only four are sufficient. It cannot be 
stated here that the elements of the structure in question are at a 
substantial distance from each other. The distance between cylinders 
is one diameter with kd^. 
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TABLE 2. TOTAL SCATTERING CROSS SECTIONS 
OBTAINED BY ITERATION AND SELF CONSISTENT 
METHODS FOR STRUCTURES CONSISTING OF TWO 
CYLINDERS: K 1 m K 2 m 90° t BfBfl 0 


xa 

xef « 4,Oxo 

xt/* r,Oxa , 

xa 

fas:" 

By 9 

Cmooomc 

BrepamoHR 

h 

xo 

Cmmcamo. 


0.2 

0.0622719 

0.0622625 

0.2 

0 . 068*6 

0.066309 

0.3 

3.439334 

3.436422 

0.3 

2.743066 

2.743112 

0.35 

1.667383 

1.667877 

0.82 

6.078781 

6.081347 

0.46 

0.661186 

0.68 IIB 7 

0.33 

4.120604 

4. 1? 1048 

0.6 

0.606346 

0.606413 

0.35 

I . 564921 

1.55185 

0.8 

0.631977 

0.532443 

0.4 

0.458834 

0.458818 

X .0 

0.551463 

0.551441 

0.46 

0.306071 

0.306863 

1. 1 

0.586818 . 

0.586090 

0.5 

0.282 M 

0.282149 

1.2 

0.643261 

0.642972 

0.7 

0.463291 

0.453373 

1.4 

0.850676 

0.860773 

0.8 

0.696771 

0.596738 

J .5 

0.995453 

0.995359 

X .0 

0.963409 

0.963266 

1.6 

1.14778 

1. 147332 

1. 1 

1.186683 

I. 186610 

1.7 

1.29306 

1.292861 

1.2 

I. 321486 

1.321451 

X .8 

1.441822 

1.441889 

1.28 

1.222493 

1.222398 


Key: a. Self consistent 

b. Iteration 

We also investigated the case when the cylinders are located 
very close to each other (but do not overlap). To calculate the 
total scattering cross section with less than 5 % error even in the 
case kd-2.1 ka, only eight iterations are sufficient (see Table 9). 
To obtain results with very high accuracy (in calculation of both 
the near and distant fields), it is sufficient to account for 18-19 
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successive scatterings. In this case, the correction to the Fourier /37 
coefficients in the 19th iteration is less than 10’"^ (see Table 10), 
which is less than 0.1/C. 

TABLE 3. HEAL PARTS OF FOURIER 
COEFFICIENTS OF SURFACE CURRENT DENSITY 
FUNCTIONS OF DIFFERENT SCATTERING ORDERS WITH 
ka-0.4, ^-52-90°, e^e^l®, kd-4.0 ka; 

PRECISE VALUES OF COEFFICIENTS: 


it*/*. - -0.0666600646. 
- 0. I360344I4, 
• 0. 243444607 


a ff 

wci) 






* fZ'/'-i 

i 

-0.3281 

-0.3251 

-0.1578 

-0.1578 

0.1578 

0.1576 

2 

0.1569 

-0.I66I 

0.0134 

-0.1444 

0.0753 

0.2331 

3 

0.0639 

-0.1024 

0.0056 

•0.1466 

0.0II6 

0.2447 

4 

0,0151 

-0.0673 

0.0010 

-0.1379 

0.0001 

0.2449 

5 

0.0017 

-0.0656 

- 0.0000 

-0.I37S 

-0.0009 

0.2440 

6 

•0.0005 

-0.0661 

- 0.0001 

0.1360 

-0.0004 

0.2436 

7 

•0.0004 

•0.0665 

- 0.0000 

-C.I360 

- 0.0001 

0.3436 

tf 

• 0.0001 

-0.0666 

- 0.0001 

-0.1300 

- 0.0000 

0.2435 


Key: a. Iteration number (i) 


Conclusion 


By using the method of solution of the problem of multiple 
scattering by a finite number of scatterers developed for systems 
which consist of smooth nonresonant components, the problem of mul- 
tiple scattering of a H polarized ws /e by N unclosed screens in the 
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form of cylinders with longitudinal slitB with distinct resonance 
properties was solved. 

TABLE IMAGINARY PORTIONS OP FOURIER 
COEFFICIENTS OF SURFACE CURRENT DENSITY 
FUNCTIONS OF DIFFERENT SCATTERING ORDERS WITH 
ka«0.H, €^2-90°, e -e 2 -l°, kd-4.0 k&; 

PRECISE VALUES OF COEFFICIENTS: 


- I.4I68I392, 

-0.151179847, 

-- 0.184379559 


p a 

m (0 

i-fi 

j-Zfr 

— n 




I 

1.3203 

1.3203 

0.1409 

0.1499 

•0.1499 

-0.1499 

2 

0.il29 

1.4532 

0.0066 

0.185’ 

-0.0154 

-0.1653 

3 

0.000b 

1.4540 

-0.0023 

.0.1533 

-0.0128 

-0.I79I 

4 

-O.0IO4 

1.423b 

•0.0015 

0.»I8 

-0.0060 

-0.B3I 

S 

-0.0060 

I.4IS6 

-0.0006 

0.J8I* 

•0.0012 

-0.1843 

6 

-0.0015 

1.4 1 VI 

-0.0001 

0.I5II 

-0.0001 

-0.1845 

7 

•0.0003 

1.4166 

-0.0000 

0.ISII 

0.0000 

-0.1844 

8 

0.0000 

1.4 MB 

0.0000 

0.1511 

0.0000 

0.1844 


Key: a. Iteration number (i) 

The problem was solved by two methods: self consistent, which 

takes full account of the interaction between cylinders, and the 
iteration. In which the degree of accounting for the interaction is 
directly proportional to the number of iterations. 

It was shown on the basis of numerical calculations that the 
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iteration method of solution gives correct results, even in the case 
of cylinders located very close to each other (more than 20 itera- 
tions must be taken into consideration here). 

TABLE 5 . REAL PARTS OP FOURIER . 

COEFFICIENTS OF SURFACE CURRENT DENSITY 
FUNCTIONS OF DIFFERENT SCATTERING ORDERS WITH 
ka=0.M, ^=^=90°, e^e^l 0 . kd=7.0 ka; 

PRECISE VALUES OF COEFFICIENTS: 


*e'/v. - - 0 . 419716396 , 

tie fa - - 0 . 168262608 , 

fit'//., - 0.196 13 1449 


as? 

Hrepa - 

IW (0 

e*jUc 



wju; 



I 

-0.3251 

-0.3251 

-0.1578 

-0.1578 

0.1578 

0.1578 

2 

- 0.0537 

- 

-0.3788 

0.0015 

-0.1564 

0.0366 

0.1944 

3 

-0.0294 

-0.4082 

- 0.0010 

-0.1574 

0.0029 

0.1973 

4 

- 0.0089 

-0.4I7I 

-0.0006 

-0.1580 

-0.0005 

O.I967 

5 

-0.0021 

-0.4192 

-0.0002 

-o. 

-0.0004' 

0.1964 

6 

-0.0004 

-0.4196 

- 0.0000 

-0.1582 

-0.0001 

0.1962 

7 

- 0.0001 

-0.4197 

- 0.0000 

-0.1582 

- 0.0000 

0.1962 

6 

- 0.0000 

-0.4197 

- 0.0000 

-0.1582 

- 0.0000 

0.1962 


Key: a. Iteration number (i) 

The basis for solution of the problem of multiple scattering by 
a finite number of smooth scatterers Is solution of the problem of 
wave multiple scattering by a single element. 

Efficient use of the iteration method became possible due to 
the solution of the problem of a single unclosed screen obtained 
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TABLE 6 . IMAGINARY PORTIONS OF FOURIER CO- 
EFFICIENTS OF SURFACE CURRENT DENSITY FUNCTIONS 
OF DIFFERENT SCATTERING ORDERS WITH 
ka= 0 .*», C 1 =C 2 *90°, e^e^-i 0 , kd= 7 . 0 ka; 

PRECISE VALUES OF COEFFICIENTS: 


. 1. 50478019, 

. 0.1677396 16, 
Z,#, • 0.0909355689 


Ap 

■wpa- 

*"«) 




$ 1 

Jmfe 


I 

1.32033 

1.0320 

0.1499 

0.1499 

-0.1499. 

-0.1499 

a 

0.I55I 

1.4754 

0.01369 

0.1636 

0.0444 

-0.1065 

3 

0.02751 

1.5030 

0.0036 

0.1670 

0.0113 

-0.0942 

4 

0.0027 

1.5056 

0.0006 

0.1677 

0.0027 

-0.0914 

5 

-0.0004 

I. 50521 

0.0001 

0.1677 

0.0006 

-0.0909 

6 

-0.0003 

r.5049 

0.0000 

0.1677 

0.0001 

-0.0909 

7 

- 0.0000 

1.5049 

- 0.0000 

0.1677 

0.0000 

-0.09068 

8 

- 0.0000 

1.5049 

- 0.0000 

0.1677 

0.0000 

-0.09068 


Key: a. Iteration number (i) 


in analytical form (in the case of narrow slits, 9<10°), i.e., just 
as for a smooth scatterer (closed screen). 

By using the procedure of inversion of the static portion of 
the matrix of the system of equations in solution of the problem of 
multiple scattering by a single cylinder with a longitudinal slit 
[21], the problem of multiple scattering by N unclosed screens by 739 

the self consistent method i*as reduced to solution of a system of 
Fredholm equations of the second kind, more precisely, to a set of 
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TABLE 7. REAL PORTIONS OP FOURIER 
COEFFICIENTS OF SURFACE CURRENT DENSITY FUNCTIONS 
OF DIFFERENT SCATTERING ORDERS WITH 
ka-1.57, ^-^* 90 °, kd=4.0 ka; 

PRECISE VALUES OF COEFFICIENTS: 


&/V. - 1,01885437, 

. -0.399932611, 
Aefa - -0. 1 1026788 


a*» 

wepa- 


Z*'fe 

kjVj 




I 

1.1883 

I. 1883 

0.0443 

0.0443 

-0.0443 

-0.0443 

2 

-0.2602 

0.9281 

-0.2333 

-0.1892 

-0.4963 

-0.5425 

3 

0.1325 

1.0606 

0.III2 

-0.0780 

0.1994 

-0.3432 

4 

-0.0585 

I.002I 

-0.0439 

-0.I2I9 

•0.0773 

-0.4204 

5 

0.0227 

1.0248 

0.0154 

-0.1066 

0.0269 

-0.3934 

n 

-0.0079 

1,0169 

-0.0047 

-0.III3 

-0.0083 

-0.4018 

7 

0.0024 

I.0I93 

0.0012 

-O.1 101 

0.0021 

-0.3097 

6 

-0.0006 

I.0I87 

-0.0002 

-0.1103 

-0.0006 

-0.4003 


Key: a. Iteration number (i) 

N simultaneous systems of linear algebraic equations of the second 
kind. 


The advantages of the iteration method of solution over the 
self consistent method were shown in the work. In our opinion, the 
most significant advantage of the iteration method is the physical 
foundations of its construction. It is evident that the degree of 
interaction between the structural elements depends on many parameters, 
primarily on the relationships between the wavelength and dimensions 
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TABLE 8. IMAGINARY PARTS OP FOURIER 
COEFFICIENTS OF SURFACE CURRENT DENSITY FUNCTIONS 
OF DIFFERENT SCATTERING ORDERS WITH 
ka*1.57, C 1 =C 2 = 90 o > e 1 «e 2 »l°, kd=4.0 ka; 

PRECISE VALUES OF COEFFICIENTS: 


Xmjtf. - 0.47I6833I8, 
r~'/U, - 1.3047184b, 
Jmfa - -1.76185287 


3 J»epa- 

r “ (0 



In ft 


JmjUj 

* 

I 

0.7863 

0.7863 

1.4847 

UA&P 

-1.4847 

-1,4047 

w*: 

2 

-0.4005 

0.3858 

-0.2180 

1.2667 

-0.3439 

-1.8286 

3 

0.1043 

0.4901 

0.0413 

1.3080 

0.0723 

-1.7563 

4 

-0.0202 

0.4699 

-0.0008 

1.3072 

-0.0012 

-1.7573 

5 

0.0000 

0.4699 

-0.0050 

1.3022 

-0.0068 

-1.7663 

• 6 

U.0027 

0.472b 

O.uOSb 

I.30b6 

0.0063 

-1.7800 

7 

-0.0019 

0.4706 

-0.0018 

1.3040 

•0.0032 

-1.7632 

8 

0.0009 

0.4717 

0.0008 

1.3047 

0.0023 

-I.76IB 


Key: a. Iteration number (i) 

cf the screens and the distance between the structural elements. For 
calculation of a multicomponent structure, it therefore is very use- 
ful to indicate zones of strong and weak interaction. The iteration 
method of solution makes it possible to do this. The iteration method 
of solution of the multiple scattering problem has still other advan- 
tages here. Because of the strict sequence of operations, it is easy 
to perform it by computer, and it does not require large memory 
volume. All this permits the solution of such complex problems even 
on a medium capacity computer. 


29 




TABLE 9- TOTAL SCATTERING CROSS SECTION 
OBTAINED BY ITERATION AND SELF CONSISTENT 
METHODS FOR STRUCTURES CONSISTING OF TWO 
CYLINDERS: 5 =go°, 0 1 =e 2 =-5° , kd-2.1 ka 


Ka 

a *4o 

CamooritacoBaHHuft mbtoa 

b 'T | 

0.2 

0.13369 

0.13102 

0.22 

0.20697 

•^T 

0.20292 

0.25 

0.40962 

m 11 W'T| 

0.40130 

0.3 

I . 5108 

1.4963 

0.31 

1.9937 

1.9904 

0.33 

3.2112 

3.2243 

0.36 

3.7107 

3.6646 

0.39 

2.4939 

2.6186 

0.4 

2.1628 

2.1325 

0.45 

1.3245 

1.3217 

0.46 

I. 1568 

1. 1549 

0.5 

1.1040 

1. 1018 

0.55 

1.0765 

1.0726 

0.6 

. I. 1221 

I. 1177 


Key: a. Self consistent method 

b. Iteration method 

In conclusion, the authors thank V.P. Shestopalov for con- 
tinual attention to the work and comprehensive discussion of the 
results . 
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TABLE 10. ZERO FOURIER COEFFICIENT OF 
SURFACE CURRENT DENSITY FUNCTION OF 
DIFFERENT SCATTERING ORDERS FOR STRUCTURE 
OF CLOSELY SPACED CYLINDERS, ka=4.0, 
e i“ e 2“ 5 °» $^ 2 = 90 °, kd-2.1 ka 


Houep a 

xrepupa 

£e/V: 

* 



I 

0.44347 

0.44347 

374469 

3.4469 

2 

0.30470 

0.74817 

-2.3135 

I. 1334 

3 

-1.0369 

-0.2888 

0.84754 

1.98094 

4 

0.75796 

0.469143 

-0.06775 

I.9I3I8 

5 

-0.38434 

0.064795 

-0.25552 

1.65766 

6 

0.06471 

0.169509 

0.24668 

1.90424 

9 

0.05248 

0.19442 

0.002346 

1.80668 

II 

0.00382 

0.17473 

0.0I74U 

1.80612 

12 

0.00444 

0.I79I7 

- 0.009636' 

1.79659 

14 

0.00352 

0.17727 

0.00060 

I. 8000 I 

15 

-0.00I4I 

0.17585 

-0.00153 

1.79847 

16 

0.000II 

0.17597 

0.00121 

1.79969 

lb 

-0.00036 

0.17595 

0.00014 

1.79923 

19 

0.00023 

0.176 19 

0.00007 

1.79930 

20 

-0.00008 

0.I76IU8 

• 0.00011 

I. 79918 


Key: a. Iteration number 
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